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Abstract 

In addition to its long-term constancy, the Pioneer (spacecraft) anomaly appears to only exist for 
bodies whose mass is less than that of: planets, moons, comets, and heavy asteroids of known mass. 
Assuming the observational evidence is reliable and not the result of an unknown systematic effect, a 
violation of the Weak Principle of Equivalence is implied. This constraint is the most confronting for 
any prospective new physics. Any new physical mechanism that proposes an additional gravitational 
force, i.e. an additional spacetime curvature, is rendered unreasonable because all masses should be 
equally affected. This paper examines an approach that is based upon the existence of a sum of ad- 
ditional field energies. A finite number of tiny wavelike undulations, upon the existing gravitational 
field, are hypothesized. The sources of these non-Einsteinian gravitational waves are the spin-orbit 
coupled moons of the solar system. An excess energy arising from: a lunar orbital motion, quantum 
mechanical geometric phase and spacetime curvature generates these new gravitational field waves. 
General Relativity's Lorentz invariance demands that these "acceleration-waves" have constant am- 
plitude. The dissipation of these spherical waves, as they expand out from suitable moons, is seen 
to exist as a volumetric based reduction in the inertial mass that can sympathetically oscillate with 
these waves. Therefore, masses above a given wave's cut-off mass remain completely unaffected. The 
full substantiation and quantification of these undulations upon the gravitational field is deferred to 
another paper. This article simply seeks to show that a superposition of these "acceleration-waves" re- 
sults in an oscillatory or non-steady expression of spacecraft kinetic energy. The superposition of these 
additional oscillatory components of longitudinal motion leads to a shortfall of actual steady trans- 
lational motion, as compared to predictions that assume a purely steady motion. With wave energy 
proportional to lunar mass, and some moons inactive for geometric reasons, the four Galilean moons 
of Jupiter and Saturn's Titan dominate this new effect. The orbital resonances of Jupiter's Galilean 
moons, markedly attenuates the variation in magnitude of the Fourier-like superposition of these waves. 
The small variance observed for the Pioneer (acceleration) anomaly around its long-term average has 
created the misleading impression of a constant additional inward acceleration with associated obser- 
vational noise. Additionally, Jupiter's least orbitally resonant moon Callisto, and Saturn's Titan, cause 
the Pioneer spacecraft to exhibit a (synodic) 356 day resonance, that has been misinterpreted as an 
Earth based 'annual' residual. The amplitudes of the cyclic diurnal and annual Earth-to-spacecraft 
motion 'offsets' are examined in some detail. The Earth based annual residual is found to be actually 
very small and incapable of figuring in the Doppler tracking observations. In order to clearly establish 
and scrutinize the physical validity of the hypothesis being presented, aspects of the author's broader 
model are incorporated into the discussion only where necessary. Direction cosine corrections, for 
wavefront direction relative to spacecraft trajectory, are quantitatively neglected in order to primarily 
scrutinize the physical viability of the hypothesis. Finally, the hypothesis sheds promising light upon 
other concerns regarding the modelling of gravitation influencing 'light' bodies in our solar system. 
These include: the "Earth Flyby Anomaly" , an apparent absence of small comets, an apparent paucity 
of smaller bodies in the Main Belt of asteroids, and residual doubts concerning the "Migrating Planets" 
hypothesis - that addresses the too rapid formation of the ice giants Uranus and Neptune. 
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1 Introduction 

The Pioneer anomaly has failed to be explained by 
any systematic effects and has two awkward obser- 
vational characteristics. The first is its (long term) 
constancy and the second is that it implies an ap- 
parent violation of the principle of equivalence — 
in that only spacecraft and not: planets, moons, 
comets or large asteroids of known mass, seem to 
be affected p. 3]. By way of these two awkward 
observational constraints the anomaly has resisted 
full explication. The latter constraint and obser- 
vational evidence appear to preclude any standard 
force based hypothesis, such as a modification of 
Einstein's General Relativity. General Relativity, 
as a theory of gravitation, holds impressively in the 
solar system for electromagnetic wave propagation 
and the motions of 'heavy' bodies. 

A hypothesis concerning the re-expression of a 
tiny proportion of total spacecraft kinetic energy 
into a number of (longitudinal) oscillatory compo- 
nents of kinetic energy is proposed. These coexist 
with and in addition to the dominant steady (non- 
oscillatory) kinetic energy. These additional 'spec- 
tral' components, together, then cause an ongoing 
shortfall in actual speed relative to predicted space- 
craft translational speed. This motion is a response 
to tiny harmonic 1 undulations of the gravitational 
field. 

1.1 Outlining a new approach 

This write up primarily presents the ideas and 
mathematics that illustrate how the anomalous 
Pioneer 'acceleration' can alternatively be seen 
as potentially a Fourier summation of first or- 
der constant amplitude (acceleration or gravita- 
tional 2 ) fluctuations upon the pre-existing gravita- 
tional field. To support this it is argued that, un- 
like the Earth based diurnal residual, the 'annual' 
residual is spacecraft based and hence real (and not 
quite of 365.25 days duration either). This residual 
is modelled to be the only obvious resonance of a 
number of sinusoidal undulations on (and of) the 

1 In order to avoid confusion with wave resonances, the 
expression 'harmonic' is not used to describe sinusoidal and 
co-sinusoidal functions. The term 'sinusoidal' is preferred. 

2 To avoid confusion with Einstein's gravitational waves, 
the expression 'acceleration waves' is preferred. These accel- 
eration waves, like gravitation, produce spacetime curvature. 



gravitational field. The amplitudes of these sinu- 
soidal undulations are very small compared to the 
strength of the gravitational field. Only in very 
weak gravitational fields is their presence signifi- 
cant. 

The source of these waves shall not be fully elu- 
cidated 3 , other than to say, that observational ev- 
idence and the author's model dictate that they 
originate from regular moons of the solar system 
held in spin-orbit coupling around their respective 
host planets (which in turn orbit the Sun) 4 . 

This new mechanism is dominated by the big- 
ger moons of the solar system. Of the big seven 
moons — Earth's moon (Moon, or unofficially 
Luna) and Neptune's Triton do not 'generate' ac- 
celeration waves upon the pre-existing gravitational 
field. Thus, Jupiter's four Galilean moons and Sat- 
urn's Titan are seen to dominate the mechanism's 
effect in our solar system 5 . 

2 Expanding the hypothesis 

In this section the basic preliminary material re- 
quired to support both: a new model, and the 
mathematics that follows, is presented. 

2.1 A Fourier wave summation 

The constancy of the Pioneer anomaly means that 
any Fourier-like summation of first order fluctua- 
tions of acceleration upon the gravitational field 
would need to be special in that: they are of con- 
stant amplitude, and together they display only a 
hint of abnormally large variance around a constant 
mean Fig. 14, p.24]. 

3 The physical establishment of the mechanism is vital, 
but this involves a long write up that is primarily quantum 
mechanical in nature. The reader is asked to 'bear with the 
hypothesis'. 

4 A fully detailed model produces amplitudes of the waves, 
and the distribution of the cut-off mass with respect to dis- 
tance from the (finite number of) sources. 

5 The collision based origin of the Earth's moon, with its 
large angular momentum, and the retrograde motion of Nep- 
tune's Triton makes their relationship between third-body 
orbital kinematics, spacetime curvature, quantum mechan- 
ical indeterminacy, and geometric phase advance incapable 
of generating any atomic virtual 'excess' energy. Since all 
atoms of some moons collectively share the same virtual ex- 
cess spin angular momentum based energy these virtual en- 
ergies 'sum' to become a singular real excess energy that is 
expressed upon the gravitational field as a new kind of wave. 
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With three of Jupiter's four Galilean moons in 
orbital resonance, and Ganymede and Callisto in a 
subtle 7 to 3 resonance; only Titan is a free 'agent'. 
These orbital resonances act to smooth out the vari- 
ations in the Fourier summation of these accelera- 
tion (or gravitational field) waves. Thus a super- 
position of waves approaches a constant amplitude 
but this is never attained in the short term. 

2.2 (Acceleration) Wave aspects 

This section is well short of fully comprehensive. 
An asymmetrical interaction of: gravitational cur- 
vature of space, the quantum mechanical geometric 
phase of atoms, and quantum indeterminacy of en- 
ergy (per orbit) in the third body of a three body 
celestial system leads to these 'acceleration' waves 
upon (and of) the gravitational field. Indeed they 
violate the usual scope of energy conservation, al- 
though, it should be noted that conservation laws 
rely upon symmetry. 

... where there is a symmetry there is a 
conservation law, and with certain reser- 
vations the converse is also true p. 159]. 

Geometric phase advance in prograde celestial 
three body motion, involving only positive mass, 
is an inherently asymmetrical situation 6 . 

General Relativity's Lorentz invariance demands 
constant amplitudes. The dissipation of the waves 
is enacted by a wave volume based reduction in 
the mass that these waves can influence to sym- 
pathetically 'oscillate' in response to them. The 
waves have: an overall spherical shape, a period 
matching the heliocentric period of orbiting moons, 
the propagation (of constant wave phase) is at the 
speed of light, and the wave's particles are neces- 
sarily spin or spin zero . The particles associated 
with these waves impart no momentum to the Pi- 
oneer spacecraft (or any other body.) Thus, it is 

6 A force based model of the union of quantum mechanics 
and general relativity always conserves energy, a priori. By 
way of an overlooked energy (imbalance and subsequent) 
transfer, e.g. from the indeterminate micro ('sub'-quantum) 
world to the macro (specetime curvature) world, this force 
based approach could be found wanting. 

7 Thus, the waves carry no angular momentum. The 
waves only 'carry' energy and yet their origin removes the 
excess virtual angular momentum per orbit shared by nu- 
merous atoms of certain moons. When such waves are gen- 
erated, the momentum of a moon, as a bulk object, remains 
unchanged. 



only the undulations upon the gravitational field 
that are capable of physically influencing ('lighter') 
bodies in the solar system. 

2.3 Regarding the equivalence prin- 
ciple 

Regarding the equivalence principle Thibault 
Damour offers the following advice. 

The Equivalence Principle (EP) is a 
heuristic hypothesis which was introduced 
by Einstein in 1907, and used by him to 
construct his theory of General Relativity. 
[The] EP is not one of the basic princi- 
ples of Nature (like, say, the Action Prin- 
ciple, or the correlated Principle of Con- 
servation of Energy). It is a "regional" 
principle, which restricts the description 
of one particular interaction [mediated by 
a massless spin-2 field]. An experimen- 
tal "violation" of the EP would not at 
all shake the foundations of physics (nor 
would it mean that Einstein's theory is ba- 
sically "wrong"). Such a violation might 
simply mean that the gravitational inter- 
action is more complex than previously 
assumed, and contains, in addition to the 
basic Einsteinian spin-2 interaction, the 
effect of another long-range field. (From 
this point of view, Einstein's theory would 
simply appear as being incomplete.) 0] 

Einstein's Principle of Equivalence is (non-locally) 
restricted to uniform fields, and thus it does not 
conflict with undulations placed upon a (predomi- 
nantly static) gravitational field that arc not solely 
of gravitational origin 8 . Indeed, a uniform gravita- 
tional field may be seen as an oscillatory field whose 
frequency goes to zero, and whose period goes to 
infinity. 

The new model being proposed, in line with the 
observational evidence, violates the weak principle 
of equivalence in that the existence or non-existence 

8 In the broadest sense, it is an interaction between gen- 
eral relativity and quantum mechanical (energy) indetermi- 
nacy (over a 'cycle' time) that is involved in the establish- 
ment of the undulations upon the gravitational field. An 
alternative name for these waves is "gravito-quantum radi- 
ation" . 
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of a body's oscillatory acceleration response to the 
field waves is (inertial) mass dependent 9 . 

2.4 Gravitation & potential energy 

In General Relativity, potential energy (P.E.) as a 
sum of particle energies is not well defined, because 
potential energy lies in the geometry of spacetime 
itself. Additionally, we may say gravitational en- 
ergy cannot be localized. The following quote from 
Michael Mensky reinforces this point of view. 

The question of conserving energy- 
momentum in General Relativity (GR) al- 
ways attracted much attention. One of 
the reasons is that covariant description 
of energy-momentum seems to be incom- 
patible with the integral conservation law. 
Particularly, it is generally believed that 
no integral conservation law follows from 
the covariant differential conservation law 
for the energy-momentum tensor (EMT) 
of matter (i.e., from its covariant diver- 
gence being zero) p. 261]. 

Today we may still think in terms of P.E. but at- 
tempting to quantitatively relate this to GR ap- 
pears to be very awkward, if not ill-conceived. On 
the contrary, relating P.E. to classical-like waves or 
'ripples' of acceleration, i.e. spacetime curvature 
undulations, is conceptually simple. The idea that 
a wave contains energy is familiar; thus an acceler- 
ation wave can be seen to possess some sort of po- 
tential energy — since it is physically a distortion 
of spacetime (in the manner of ripples on a pond) . 
The existence of these field undulations makes use 
of the curved spacetime conceptualization of GR, 
but at non-special-relativistic speeds. 

From the above, it appears that the understand- 
ing of the links between GR and energy (and grav- 
itation) may be accepted as incomplete. 

2.5 Gravitational field oscillations 
and general relativity 

Sections 12.31 and 12.41 sought to show that there 
appears to be nothing in the foundations of gen- 

9 Additionally, the wave's energy is dependent upon the 
number of atoms or molecules in a moon and hence the com- 
position of a moon shall also have a slight effect upon the 
total amount of wave energy a moon (of given mass) may 
'release' or 'generate'. 



eral relativity, itself, that distinctly forbids a new 
(gravito-quantum) mechanism from generating os- 
cillatory fluctuations upon a pre-existing gravita- 
tional field 10 . When a moon rotates around a 
planet a similar oscillatory field effect, although 
non-constant with radial distance, is experienced 
by a point mass — in the not too near vicinity. 

The generation of the waves that affect Pioneer 
spacecraft (S/C) are quite distinct from three body 
effects and general relativity's gravitational waves. 
They may be thought of as type-2 gravitational 
waves if one prefers. Indeed, the existence of type-1 
gravitational waves indirectly supports the exis- 
tence of another kind of gravitational wave 11 . 

2.6 Cut-off mass, convolution, and 
wave energy dispersion 

A sharp (all or nothing) cut-off 12 exists for the ef- 
fect of an individual acceleration wave upon masses 
in motion. 

Inertia involves a body's ability to resist a change 
in motion. With the hypothesis being presented, it 
is also necessary to see inertia as associated with a 
wave field's ability to 'instill' an oscillatory varia- 
tion in the motion of a body. 

This situation may be aligned to the mathemat- 
ical technique of convolution, which determines a 
system's output [oscillatory behaviour: yes or no] 
given an input signal [the acceleration wave] and 
the system impulse response [some function of the 
inertial mass of a body] 13 . 

Although (acceleration) wave amplitudes remain 
constant, the (inertial) mass associated with this 
new kind of wave reduces in proportion to the vol- 
ume the dispersing wave encompasses 14 . Thereby, 
once a type-2 gravitational wave is generated (or 
established) , conservation of energy is obeyed 15 . 

10 Certainly, the principle of special relativity is in need 
of some attention, and the spirit of GR seems to be under 
threat, but the provisional acceptance of these type-2 grav- 
itational waves is not unreasonable. 

11 A further source of non-uniformity for a gravitational 
field is tidal effects. Such effects ensure lunar spin-orbit 
coupling around the moon's host planet. 

12 Much in the manner of the Photoelectric Effect. 

13 Paraphrasing: http://www.see.ed.ac.Uk/~ mjj/dspDemos 
/EE4/tutConv.html (square brackets content excluded). 

14 Once again deference to a fuller elucidation of the model 
is necessary. 

15 It is only the generation of the acceleration- waves that 
appears to disobey (expected) conservation of energy. 
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3 Amplitudes associated with 
acceleration waves on the 
gravitational field 

This section begins the process of quantitative sup- 
port for the hypothesis. Two results are estab- 
lished. Results pertaining to the relationships be- 
tween the amplitudes of sinusoidal: acceleration, 
speed and range variations of a body are outlined. 
Since the wave amplitudes (via Lorentz invariance) 
are necessarily constant, the mathematics that fol- 
lows is greatly simplified. Note that the solar sys- 
tem barycenter acts as the (quasi-) 'global' (i.e. so- 
lar system) inertial frame's reference point. 

3.1 Amplitudes of the cyclic varia- 
tions 

The integral of a constant acceleration (of unit am- 
plitude) over a time ^ is simply a velocity of magni- 
tude ^ . Now J Q 2 sinO d9 — 1 is a velocity amplitude 
related to a unit amplitude sinusoidal acceleration 
acting over -| (i.e. a quarter of a wavelength). The 
period of the wavelength may be either 2ir or At. 
Closely related to the maximum sinusoidal acceler- 
ation (wave) amplitude Aa and a time of At/4 is 
the maximum amplitude of a similarly sinusoidal 
velocity so that: 
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Notice that we let: 

| A a*! = Aa, |Alf| = Av and |Ax*| = Ax. 

In short, there is a direction and time indepen- 
dent relationship between the magnitudes of veloc- 
ity and acceleration wave amplitudes 16 . This ap- 
pears in equation 50 of Anderson et.al. |2J p. 37] 
with Aq replacing Aa. 

Similarly, we also have for sinusoidal speed varia- 
tion and associated (j wave length) sinusoidal range 
change, the following relationship of amplitudes: 
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16 By treating magnitudes only we are effectively now only 
discussing spacecraft motions beyond Saturn, with trajec- 
tories assumed orthogonal to the wavefronts. Effects from 
the smaller moons orbiting Uranus, Neptune and Pluto are 
quantitatively negligible. 



Note that: Wdiumai ~ 7.3 x 1CP 5 rad/s and 
^annual ~ 2.0 x 10 -7 rad/s. Additionally, by way 
of Ref. [2, pp.8, 15, 37], the (return trip) Doppler 
frequency shift {Av) is determined via: 



Av 
v 



2 dl 
c dt 



with the S-band downlink frequency (v) being 
-2.29 GHz. Finally, for S-band Doppler: 1 Hz 
corresponds to 65 mm/s P-18], or more prag- 
matically, 10 mHz corresponds to 0.65 mm/s. 

3.2 Further comments 

Observe that if undulations in the gravitational 
field are the cause of fluctuations in spacecraft 
(S/C) speed then (assuming pure radial motion 
for the S/C and a position well beyond Saturn): 
A~ct — —LuAlf and Alt = cuAlt. 

Note that direction cosines are rampant in this 
new approach involving spherical waves emanating 
from a moon. Small corrections in the form of di- 
rection cosines, for the trajectories of bodies nearer 
the center of the solar system, are (quantitatively) 
neglected in this write up. 

The overall write up is both quantitatively ide- 
alized and theoretically lacking full substantiation, 
in order to comprehensively establish the physical 
validity of the acceleration-wave hypothesis. 

3.3 A set of annual reference values 

The results of Section 13.11 applied to an annual 
residual allow a set of reference values to be estab- 
lished. For Av — 0.2 mm/s at Av 3.1 mHz (a Pi- 
oneer S/C S-band Doppler frequency change), the 
approximate values of range and acceleration am- 
plitude are: Ax — 1 km and Aa = 0.4 x 10 _8 cm/s 2 
respectively. These values are all linearly scalable 
for different magnitudes of the four physical 'quan- 
tities' involved (Ax, Av, Av, andAa). 

4 Interpreting the Pioneer 10 
diurnal residual 

This section shall begin to examine the claim of 
Anderson et. al. that: "[the] annual and diurnal 
terms are very likely different manifestations of the 
same modelling problem 2, p. 38]." 
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4.1 Introductory remarks 

There are three causes of residuals arising from a 
method of least squares analysis. These are: obser- 
vational error, approximation of parameters, and 
model or theory inadequacy [HI Ch.8]. Doppler data 
cannot clearly discern whether an oscillatory resid- 
ual is Earth or spacecraft based. Only the relative 
acceleration, motion, or (line segment) distance be- 
tween the two bodies is determined by the mea- 
surements. Assuming the observations are reliable, 
either an oscillatory sinusoidal residual is a result 
of Earth based parameter error(s), or failing that, 
it is an unlikely real spacecraft motion and hence 
beyond current gravitational theorization. 

4.2 The magnitude of the diurnal 
residual 

The diurnal residual's interpretation is crucial to 
understanding the annual residual. The noise of the 
diurnal residual is greatest around solar conjunc- 
tion, but at solar opposition, near a minimum in 
the solar cycle 17 , exceptionally good data is avail- 
able: see Fig.18, p.38]. 

Markwardt gives a figure of lOmHz (i.e. 
0.65mm/s) for the amplitude (on average) of the 
diurnal residuals. Ref. Fig.18, p.38] implies 
a Nov/Dec 1996 solar opposition amplitude of 
0.1376 mm/s that may be obtained from: 

Aa = luAv 

and the values of uid.t. — 7.2722 x 10~ 5 rad/s and 
a d , t , = (100.1 ± 7.9) x 10- 10 m/s given by Ref. 
p.38]. This gives, via Aw = loAx, a cyclic diurnal 
position offset of ~ ±1.9 metres (i.e. Aa; « 1.9 m). 
As an angular offset at the Earth's surface this 
equates to: 

1 89 

= tan" 1 637gx — 3 « 1.70 x lO^deg. « 60mas 

where 'mas' is milliarcseconds. This being a small 
and non-problematic distance or angle, that is less 
than 0.1% of maximum DE405 ephemeris error in 
the Earth's orbit of 2 kilometres 0, and about 
16 times 1997 Earth orientation polar motion (root- 
mean-square) calibration accuracy of ~12 cm 0. 

17 There was a broad minimum in the solar cycle around 
May 1996. 



4.3 Diurnal residual parameter 
groups 

Three main groups of factors or parameters affect 
the diurnal residual. For a diurnal effect only short 
term parameters and location aspects are included 
in the first two groups. 

• Earth orientation parameters (EOP) errors. 
This concerns celestial polar motion offsets, 
and variability of the earth's rate of spin. Let 
us also include here: antenna location errors 
for NASA's Deep Space Network (DSN). 

• Planetary ephemeris errors. 

These concern the location of the earth in its 
orbit (relative to the solar system barycenter). 

• Error via miscellaneous effects. 

Including: ocean tides, weather, and vari- 
able atmosphere effects; troposphere and iono- 
sphere effects (causing spectral broadening of 
the carrier wave frequency); and interplane- 
tary scintillation (i.e. plasma-based fluctua- 
tions). 

All of this information, and much more, goes 
into either JPL's Orbital Determination Pro- 
gram (ODP) or the Aerospace Corporation's 
Compact High Accuracy Satellite Motion Pro- 
gram (CHASMP). The presence of errors or inaccu- 
racies results in a residual, or what may be termed 
a "spacecraft motion offset" . Anderson et. al. have 
found that individually the first two sources of er- 
ror cannot be (solely) responsible for the residual 
p.36]. 

The diurnal residual of Fig.18] may be said 
to contain two aspects. Firstly, a stochastic (or 
random) aspect and secondly, a cyclic aspect. In 
the author's opinion, the miscellaneous effects (dis- 
cussed above), and Earth spin rate variability, pre- 
dominantly produce either: very small effects or 
random shifts in the residual, and hence they may 
be neglected from an account of the cyclic residual 
indicated in Fig.18. It is the cyclic diurnal signa- 
ture that shall now concern us, as it most closely 
relates to the annual residual's interpretation. 

Note that the period of the diurnal cycle is not 
a concern. A diurnal residual should exist in the 
Doppler data. The Earth is a "wobbly platform" . 
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4.4 An interpretation of the cyclic 
diurnal residual's amplitude 

It appears that only a combination of EOP and 
Ephemeris error together, in the ODP or CHASMP, 
can (primarily) produce the cyclic diurnal ampli- 
tude error observed at opposition 18 . 

This interpretation is analogous to how celestial 
pole offsets arise. Celestial pole offsets are required 
because the model (of the Earth's orientation) that 
combines precession and nutation relies on fixed pa- 
rameters for the Earth's shape (geodesy) and inter- 
nal structure, but since these are not fixed the off- 
sets necessarily arise. Similarly, combining EOPs 
and the ephemeris in an orbital determination pro- 
gram is seen to produce a (pure) cyclic diurnal 
residual. This being at 0.1% of the level of the max- 
imum error in the DE405 ephemeris (i.e. 2 km). 

This seems to indirectly agree with errors in 
EOPs changing the value of a p only in the 4th 
digit p. 36]. The diurnal residuals indicate EOPs 
changing the S/C position location (only) in the 4th 
digit of the (Earth position) ephemeris error. 

4.5 In summary 

It appears that a combination of many parameters 
produces the diurnal residual, but a combination 
of two parameter 'groups' (EOP and ephemeris er- 
rors) dominates the production of the cyclic aspect 
of the diurnal residual — about a mean value. 

5 Examining the Pioneer 10 
annual residual (1987-1998) 

An understanding of the diurnal residual allows us 
to now closely examine the cyclic 'annual' residual. 
The amplitude of the diurnal cyclic residual was 
shown to be primarily a combination of errors in 
two parameter groups. Anderson et. al. claim the 
annual residual is probably due to a similar com- 
bination of parameter errors. These being: Earth 

18 This conclusion is based upon email correspondence 
with E. Myles Standish of JPL regarding the diurnal resid- 
ual. Note that Myles, whilst having no objection to this 
interpretation of the cyclic residual, emphasized the need to 
not overlook the other effects mentioned. Also see footnote 
125 of Ref. p.49]. 



based location and (long-term) polar orientation er- 
rors, and errors in the navigation program's 'deter- 
mination' of spacecraft (S/C) orbital inclination to 
the reference frame being employed. 

Note that when averages over a full day or a num- 
ber of days are used, the diurnal (and short-term) 
parameter errors, disappear from the data. Thus, 
the annual residual's amplitude is independent of 
the amplitude of the cyclic diurnal residual. 

5.1 Orbital inclination errors and 
the annual residual 

Jet Propulsion Laboratories (JPL), The Aerospace 
Corporation, and Craig Markwardt all used DE405 
in their 'best' analyses. Interestingly, beginning 
with DE400 (development ephemeris 400), both 
the Earth orientation parameters and the Earth's 
ephemeris are aligned to the ICRF (International 
Celestial Reference Frame). Error in the orbital in- 
clination of a spacecraft, via an orbital determina- 
tion program, is thus more closely related to EOP 
error than previously. 

A combined outer solar system spacecraft incli- 
nation angle error with Earth polar orientation an- 
gle error, is seen to interact with Earth location 
(ephemeris) error to produce the annual residual. 
This account of the cyclic annual residual, like the 
cyclic diurnal residual, is based on two primary pa- 
rameter groups. This account although not specifi- 
cally stated, is alluded to in Ref. p. 23, pp.36-38]. 
It appears a reasonable, although sketchy, explana- 
tion. This scenario is now further investigated. 

5.2 Quantifying the error associated 
with orbital inclination 

Standish 8, p. 1166] by way of the recent (2003) 
ephemeris DE409 has been able to quantify er- 
rors in earlier ephemerides. For all the outer plan- 
ets (short-term) inclination errors for both geocen- 
tric right ascension and declination are . 1 (arc- 
seconds) for DE200, and 0".05 for DE405. 

Since the Orbital Determination Programs 
(ODPs) 'produce' the inclination error, it is best, 
in this case, to use the error from DE200 which 
covered 1979 to 1997. In the case of a pre- 
dominantly radial motion based variation in the 
Doppler, the effects of S/C inclination error may 
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be determined 19 . Assuming radial motion for the 
distant Pioneer S/C, a (plus or minus) 0.1 = in- 
clination error for a spacecraft at 55 AU (Astro- 
nomical Unit) implies an uncertainty in position, 
orthogonal to the spacecraft's trajectory, of: 

z = (55) (150 x 10 6 )tan(0".l) w 4,000 km 

noting that tan(0".l) « 4.85 x 10 -7 . This orthog- 
onal uncertainty may then be related to a line of 
sight uncertainty by a similar triangle, such that: 

Ax = 4 x 10 6 tan(0T) « 2 m 

where Ax is the magnitude of the range variation. 
Thus, for 0".l: z w 4,000km and Ax w 2 m; and 
for l": z 40,000 km and Ax » 0.2 km or 200 m, 
whereas for 5 " : z » 200, 000 km and Aa; « 5 km. 

Thus, on their own the S/C inclination errors 
arising from DE 200 fail to account for the 'annual' 
residual, by a long way. At its minimum the 'an- 
nual' range variation amplitude (As) of Pioneer 10 
is approximately 500 meters via Av « O.lmm/s 
(see Section I3~3l and Ref. p.38]). Total orienta- 
tion errors in excess of 1.5 arc-second appear to be 
required, i.e. above 15 times the reported orienta- 
tion error of DE 200, and 30 times greater than the 
DE405 orientation errors of the outer planets. 

5.3 Remarks on combining orbital 
orientation and ephemeris errors 

It should be noted that errors in the Earth's or- 
bital position (including heliocentric radius error) 
will have a minimal cyclic impact upon the errors 
discussed in Section l5~2l Since the Earth and Pi- 
oneer 10 both lie very near the plane of the eclip- 
tic, only errors in the Earth's orbital position far 
away from Sun-Earth-Pioneer 10 conjunction or op- 
position will yield a (line-of-sight) Doppler range 
residual. This residual would show narrow peaked 
maximum or minimum amplitudes, rather than the 
smooth 'sinusoidal' wave observed |131 Fig. IB]. 

Finally, with the earlier DE 200, the ephemerides 
were oriented onto their own inherent Earth's mean 
equator and dynamical equinox. Thus, the addi- 
tional orientation error related to the Earth's an- 
nual orbital motion may slightly increase the error, 

19 This, and what follows, is related to the discussion on 
p. 74 of an article by W. G. Melbourne 1 101 on "space naviga- 
tion" , where range error is related to an orthogonal distance 
(and hence angular) error. 



but this increase would need to be larger than the 
error in the outer planets and this is unrealistic. 

Subsequently evidence to the contrary of the 
Ref. [2] stance on the annual residual is worth con- 
sidering (this is pursued in Section EJ. Addition- 
ally, let us note that Markwardt [7J p. 11] believed: 
"... the source [of the annual residual] was ulti- 
mately undetermined." 

5.4 A stance on the 'annual' cyclic 
residual 

The annual spacecraft motion offset, and hence the 
position offset, obtained from Doppler tracking ob- 
servations — appears to be neither: a purely orbital 
inclination effect, nor is it feasibly a combination of 
this with either: Earth orbit orientation error, or 
Earth position error. 

The reason for the 'annual' residual appears to be 
restricted to a choice between: a mixture of Earth 
position (ephemeris), Earth orientation, and space- 
craft inclination errors; or alternatively, an inade- 
quate model of the spacecraft's motion. Concerns 
regarding the validity of an account based on pa- 
rameter errors were raised in Sections 15.21 and 15.31 
Indeed, the "orbital determination programs" need 
to be of a high quality (in the first place) to unam- 
biguously obtain the Pioneer anomaly. 

6 Reinterpreting the Pioneer 
spacecraft 'annual' residuals 

Section found that spacecraft orbital inclination 
error appears insufficient on its own, or together 
with Earth ephemeris and orientation errors, to ac- 
count for the 'annual' residual's amplitude mea- 
sured by Doppler tracking observations. Subse- 
quently, the residual appears to be due to either: an 
unrealized ephemeris error of the Earth's position 
in its orbit, or the spacecraft has an unmodelled 
(i.e. real) annual longitudinal oscillatory motion. 

Further concerns may be raised regarding firstly, 
the cyclic amplitude of an Earth based explanation 
of the 'annual' residual of the anomalous Pioneer 
'acceleration'; and secondly a new concern regard- 
ing the period of this residual is raised. Beginning 
with the latter, these concerns are now addressed, 
and linkages to the acceleration-wave hypothesis 
are drawn upon and incorporated in the discussion. 
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6.1 An alternative approach to the 
Pioneer 'annual' residuals 

On page 38 of the extraordinarily comprehensive 
Physical Review D paper [2] discussing the Pio- 
neer anomaly, the angular velocity is given (in in- 
terval III) as ui a .t. = (0.0177 ± 0.0001) rad/day 
which equates to 355 ± 2 days per (2tt) cycle. Fig 
2. of Scherer et. al. jTJ shows the real part of 
the autocorrelation function of the later Pioneer 10 
data (1987-1995). By averaging, from the graph 
the clear maximum and minimum range of values, 
at the half and full year, a period of ~ 355 days 
(and not ~ 365 days) is confirmed. The shape of 
the real part of the autocorrelation function indi- 
cates a solitary sinusoidal-like oscillation dominates 
the spectral aspect of a time series representing the 
(long-term) Pioneer Doppler data. Markwardt 13 
p. 11] refers to the "^annual" residuals, also identi- 
fying a mean period that is far enough away from 
365.25 days (i.e. 3%) to be worthy of signification. 

The (heliocentric-based) orbital periods of 
Jupiter's moon Callisto and Saturn's Titan are re- 
spectively: 16.689018 and 15.945421 days. Re- 
markably, their periods will resonate every 357.9 
days (m = n + 1 where n w 21.445) 20 . By making 
a synodic correction for the location of the Pioneer 
10 spacecraft with respect to the motions of the 
host planets of these moons (1992.5 to 1998.5), the 
period goes to approximately 356.1 days 21 . This 
cycle has a significant resonance amplitude, and is 
the only one freely visible, involving the (proposed) 
lunar generated acceleration-waves. Remember, if 
output from only five moons dominates the Pioneer 
anomaly, then only the Callisto-Titan resonance 
is expected to be significantly under-affected by 
the orbital resonances of Jupiter's Galilean moons. 
These resonances act to minimize the (statistical) 
variance of a p observations through time 22 . 

20 For planetary ring systems at least: "Resonances are 
strongest when m = n + 1 (for example 2:1 or 43:42) and 
weaken rapidly as m and n differ more and more [El p.70]." 

21 In the 6 years of interval III (1992.5 to 1998.5) Jupiter 
tracks, with respect to Pioneer 10's location and trajectory, 
approx. +11° prograde, whereas Saturn's position remains 
essentially unchanged. Callisto's prograde progression is 
thus ~ 1.8° (relative to Titan) per 357.9 day resonant cy- 
cle. This yields a shortening of the (360°) resonance cycle 
of ~ 1.8 days. (See Section 16.51 for the data source used to 
establish these angles.) 

22 A small amplitude, directional cosine based, 11.86 year 
Jupiter cycle should be present, and possibly evident, in the 



6.2 Pioneer data and the ephemeris 

Any inconsistency that exists between the Pi- 
oneer anomaly's annual residual amplitude and 
DE405 ephemeris error in the earth's orbital loca- 
tion of ±1-2 km [HI p. 1171] is removed if the 'an- 
nual' residual is deemed real. Markwardt [7| ob- 
serves a lOmHz annual amplitude 23 which implies 
a 3.25 km range error (see section 13.31 for reference 
values). Turyshev et. al. pi] quote an amplitude of 
1.6 x 10 _8 cm/s 2 implying a range error of 4 km 24 . 
Only later in Interval III can the annual residual's 
magnitude be said to be 'within expectations'. 

Similarly, the month of data given by Anderson 
et. al. for the diurnal residual Fig. 18, p. 38] 
covers about 30° of the annual anomaly's cyclic pe- 
riod. An amplitude increase of about O.lmm/s (to 
a maximum) over the 30 days is evident. If this 
implies [via (1 - COS30 )" 1 = (0.134)- 1 » 7.5] a 
speed sinusoid amplitude of ~ 0.75 mm/s then this 
roughly agrees with Markwardt 's value of lOmHz or 
0.65mm/s — for the amplitude of this prospective 
two- wave resonance 25 . 

6.3 Pioneer data and pulsar timing 
experiments 

Chandler |141 p. 108] states that: the annual sig- 
nature of the Earth's motion dominates the varia- 
tions of (long term) pulse arrival times. Addition- 
ally, pulsar timing measurements are accurate to 
about 3 x 10 -6 seconds or 90 meters, and there has 
been no indication from pulsar timing experiments 
of any overly large error in the Earth's ephemeris 
(i.e. orbital location over time) e.g. pp. 718-19]. 

The non-problematic account of a slight change 
in a pulsar's location by changing from DE200 to 
DE405 where only a 0.2 mas mismatch is ap- 

long term observational data of the anomalous acceleration 
[%>(£)]. This is perhaps partially evident in |2] Fig. 14, p. 24]. 

23 Markwardt finds the rms residuals of all the 'non- 
extreme' Doppler data to be of order 8mHz — see his Ta- 
ble II. The annual residuals are thus of the same order of 
magnitude as the noise, although the noise amplitude, de- 
pending upon space and atmospheric Doppler transmission 
conditions, is quite variable. 

24 The first oscillation in 1987 has an amplitude of 2.5 X 
10 — 8 cm/s 2 indicative of a ± 6.25 km range variation. 

25 This is awkward because it is inconsistent with Ref. 2 's 
stated amplitude (0.1053 mm/s) for the interval III annual 
sinusoid, or 1.5mHz implying a ± 0.5 km range sinusoid. See 
Section 16.61 for further discussion. 
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parent, further indicates an absence of cphcmcris 
accuracy concerns. A similar comment is made by 
Anderson et. al. regarding planetary (and space- 
craft) ephemeris error p. 36]. 

Note that pulsar timing data lack the accuracy of 
Doppler diurnal residual amplitude measurements, 
that are approximately 2 meters (or ~ 0.14 mm/s) 
around the 1996 Pioneer-Earth-Sun (solar) opposi- 
tion. Also note that Doppler data precision reduces 
as the period of oscillations becomes longer (given 
a fixed oscillation amplitude of frequency variation 
Aj/). 

6.4 The Pioneer 10 and 11 'annual' 
residual amplitudes 

A combination of primarily spacecraft inclination, 
and also Earth orientation and ephemeris errors has 
been proposed as a reason for Pioneer ll's greater 
'annual' residual amplitude c.f. Pioneer 10 
p. 37]. Alternatively, this feature may be related 
to the inclination of the spacecraft trajectories rel- 
ative to the equatorial planes of Jupiter and Sat- 
urn. The planet's equatorial plane is essentially the 
plane within which Jupiter's four Galilean moons 
and Saturn's Titan orbit. Jupiter's and Saturn's 
equatorial planes are tilted at 3.1 and 26.7 degrees 
respectively, relative to their orbital planes, with 
their orbital planes inclined at 1.3 and 2.5 degrees 
respectively, relative to the plane of the ecliptic. 

Disregarding the inclinations of the planets' or- 
bits, the maximum possible Callisto-Titan reso- 
nance amplitude, from lunar generated acceleration 
waves, will thus be approximately at 15 degrees. 
The Pioneer 10 and 11 spacecraft trajectories are 
inclined at approximately 3 and 16 degrees respec- 
tively, to the plane of the ecliptic. Thus, Pioneer 
11 would be expected to have a greater 'annual' 
cyclic residual amplitude, and probably a slightly 
different period to that measured by Pioneer 10 — 
if the hypothesis being outlined in this write up is 
viable. 

If the lunar based acceleration- wave hypothesis is 
viable then the other (primary) orbital resonances 
of Jupiter's moons will probably play some role in 
varying the amplitude of this (Callisto-Titan) res- 
onance over time. This arises because the other or- 
bital resonances 'complicate' the simple decompo- 
sition of this two- wave resonance out of the overall 
superposition effect of all the acceleration waves. 



6.5 The fine details of the 'reliable' 
Pioneer observations 

A number of details regarding the Pioneer space- 
craft observations appear to be illuminated by the 
"acceleration- wave hypothesis". All that is re- 
quired is to assume the time-averaged observations 
are reliable and slightly more accurate than gener- 
ally considered. Wave motion to spacecraft motion 
direction cosines are implied in all of the following 
fine detail aspects. 

1. Why the overall values of Pioneer 10 and 11 
differ slightly 26 . 

2. Why the 'annual' residual's amplitude is 
greater for Pioneer 11 c.f. Pioneer 10. 

3. Why the maximum anomalous acceleration for 
Pioneer 10 is in early 1998. 

4. Why the Pioneer 11 anomaly increases rapidly 
post Saturn encounter [3 Fig. 7, p. 18]. 

5. Why the magnitude of the Pioneer 10 anomaly 
is slightly greater in the later data, as com- 
pared to earlier data. (See Section IV. 121) . 

6. (Possibly) Why Pioneer 10 spin rate decreases 
whereas Pioneer ll's increases between ma- 
noeuvres 27 . 

Point 4 is also influenced by wave direction. At 
any point within Saturn's orbit, waves may arrive 
from Jupiter's Galilean moons and Saturn's Titan 
at obtuse angles to each other. Waves from op- 
posing directions will act to cancel each other out 
somewhat. (Section [7| shall clarify this assertion.) 

Regarding point 3: from the Pioneer 10 space- 
craft's perspective, in early 1998 Jupiter and Sat- 
urn are closer in the sky to its '(reverse) trajectory 
line' than at any other time. They lie near the 
spacecraft's reverse (or negative) trajectory line 28 

26 Noting that Pioneer 10 has over 11 years of high quality 
data, whereas Pioneer ll's data, with only 3| years, has had 
insufficient time to establish a representative 'longer-term' 
average. 

27 Naturally gas leaks are a likely cause but "for the Pio- 
neers there were anomalous spin-rate changes that could be 
correlated with changes of the exact values of the short term 
a p . The correlations between the spin-rate changes and a p 
a were good to 0.2 X 10~ 8 cm/s 2 and better El P-4019]." 

28 Essentially this is the view (for a forward looking 'pas- 
senger') in an imaginary rear view mirror on the spacecraft. 
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and approach a conjunction. Hence direction co- 
sine adjustments are minimized. Saturn actually 
crosses the extended (negative) trajectory vector 
at, or very near to, the early 1998 maximum. See 

P-24, Fig. 14] as well as planetary and spacecraft 
positions given by way of the "National Space Sci- 
ence Data Center" website — to obtain the solar 
ecliptic reference frame coordinates through time 29 . 

6.6 Damping of the annual residual 

Fig. IB in Ref. shows the annual residual, as 
does Fig. 13 in Ref. [1 p.23]. The latter figure is 
divided into intervals: I, II, and III. Thus, it is rea- 
sonable to also break Fig. IB into three intervals. 
The annual anomaly is most impressive in intervals 

1 and II (Jan 1987 to July 1992), with a single si- 
nusoidal line able to smoothly connect essentially 
every point (not quite the line shown). Such is 
not the case with Interval III (July 1992 to July 
1998) which also contains amplitudes that are de- 
cidedly smaller than in intervals I and II. Thus, 
a fairly clear distinction exists either side of July 
1992. The interval III data is also decidedly less 
sinusoidal than interval I and II data. Unfortu- 
nately, the qualitative detail of Pioneer ll's annual 
anomaly is not discussed in Ref. 0- 

A pure two- wave resonance based sinusoidal vari- 
ation should not exhibit monotonic damping in its 
[a p (t)] time series, but there are a number of mit- 
igating circumstances. Firstly, the orbital reso- 
nances of Callisto with Jupiter's other 3 Galilean 
moons, which are not rigidly fixed resonances 30 , 
will be influential upon the Callisto-Titan reso- 
nance. Particularly prominent in the damping 
could be the method of data smoothing employed. 
There may also be another reason, unknown to 
the author, for this apparent discontinuity in data 
'quality' between intervals II and III. Curiously, the 
data in interval III is actually preferred to I, II and 

29 See http: / / nssdc.gsfc.nasa.gov/space/helios/planet.html 
for planet coordinates and .../helios/heli.html for Pioneer 
spacecraft coordinates. An enlargement of Fig. 3 in Ref. |2] 
p. 5] provides a useful adjunct to this data. 

30 Even though Io, Europa and Ganymede are implied 
as being in an exact 4:2:1 resonance, their orbital peri- 
ods lag/lead relative to each other over time. This full 
lead/lag 3 moon cycle takes on average about 1.2 years (i.e. 
~ 14.4 months.) This duration is based upon "Sky and Tele- 
scope" magazine's pictorial presentation of Jupiter's moons 
over twelve years (1987 to 1998). 



Pioneer 11 data by Anderson et. al. |21 p. 26] be- 
cause the match between JPL's ODP / Sigma and 
The Aerospace Corporation's CHASMP data is so 
good. Finally, it is worth mentioning that it is hard 
to understand how an annual residual based upon 
orbital inclination (angle) error could strongly di- 
minish with increasing S/C distance from the Earth 
- even taking into account a reduction in Earth 
and spacecraft orientation errors over time. 

Nevertheless, even with this rationalization the 
damping is a concern for the hypothesis being pre- 
sented in this write up. Particularly in interval III, 
and to a much lesser degree in intervals I and II. 
This situation could be clarified by a re-processing 
of earlier pre-1987 data [IB] [T3|. 

6.7 'Annual' residual — final remarks 

Bearing in mind Section 16 . 61 there is still an assort- 
ment of evidence that the 'annual' residual can- 
not be a result of parameter error (s), as the diur- 
nal residual surely is. Subsequently, the ^annual 
residual appears indicative of a modelling oversight. 
This model inadequacy suits the new hypothesis 
presented, with a Callisto-Titan acceleration wave 
'resonance' causing a real oscillatory motion of the 
S/C relative to an 'accurately' positioned Earth — 
itself being relative to the solar system barycen- 
ter. The single sinusoid-like variation in the au- 
tocorrelation data of a p |11| . implies that only an 
(~ 355 ± 2) day residual is measured (in Inter- 
val III spanning 1992.5 to 1998.5), implying the 
(365^ day) 'true' annual residual is too small to be 
detected. The small changes between DE200 and 
DE405 for both: the Pioneer anomaly, and pulsar 
timing experiments, add support to this assertion. 

It is unlikely a spectral analysis of a very long 
a p time-series would show any sign of the lesser 
'true' cyclic annual spectral component quantified 
for inclination error in Section f5.2l at approximately 
2 meters. Following the diurnal case fSection l4.2f) . 
let us estimate total annual error at < 40 meters. 
This signature absence is due to Doppler data pre- 
cision diminishing as the period of oscillations be- 
comes longer (given a fixed oscillation amplitude 
of frequency variation Av). For the yearly oscilla- 
tions, a O.Olmm/s amplitude, which is around the 
Doppler tracking's best level of accuracy, represents 
a 50 meter range amplitude (recalling Section EOf . 

The greater precision of a 21st century mis- 
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sion, specifically designed to test the Pioneer 
anomaly |T^], would be expected to fully clarify this 
issue. The inclusion of range measurements would 
provide a cross-reference for closely examining, and 
unambiguously explaining, this solitary ~ annual 
(S/C based), or 'true' annual (Earth based) cyclic 
residual present in the Pioneer 10 and 11 data. 

7 Gravitational field undula- 
tions and energy transfer 

In light of sufficient concerns regarding the cur- 
rent account of the annual residual of the Pioneer 
anomaly, we return to pursuing the (acceleration- 
wave) hypothesis sketched earlier. This section 
seeks to show that if the (355 ± 2) day cycle is real 
and spacecraft based (at ~ 356 days) — the origin 
and explanation of the Pioneer anomaly is open to 
a promising and progressive alternative interpreta- 
tion. 

7.1 Introduction 

It has been mentioned by Anderson et. al. p. 39] 
that a sinusoidal speed variation cannot contribute 
to the Pioneer spacecraft anomaly. This is true for 
an Earth based and/or spacecraft inclination based 
residual, e.g. by way of ephemeris errors and Earth 
orientation parameters (EOPs) errors in the orbital 
determination program (ODP); but this stance may 
be challenged if the oscillatory motion is real — and 
applied to 'light' bodies with non-zero mass 31 e.g. 
the Pioneer spacecraft. We entertain and examine 
this possibility by way of the interaction between 
spacecraft 'geodesic' motion, and the kinetic and 
potential energies involved. Recall, a geodesic is 
the closest thing there is to a straight line in curved 
spacetime. When (type-2) undulations exist on the 
gravitational field, circumstances are different from 
those of a 'static' gravitational field, especially as 
regards geodesic motion. 

With no local restoring 'force ' to oppose the ef- 
fects of a cyclic gravitational undulation, light bod- 
ies in celestial geodesic motion 'go with' the (curved 
spacetime based) cyclic undulations; whereas the 

31 Note that the overall propagation speed of electromag- 
netic radiation is necessarily not retarded. It ensues that 
without inertial mass (by definition) the kinetic energy of 
E/M radiation particles cannot be incrementally 'eroded'. 



greater inertia of 'heavy' bodies means they fail to 
respond. For light bodies, spacecraft speed and dis- 
tance will oscillate around equilibrium values. A 
standard barycenter systemic 32 inertial frame re- 
mains valid. With the passing 'phase' of an accel- 
eration wave at the spacecraft, the local gravita- 
tional acceleration field strength varies cyclically, 
relative to the system's barycenter inertial frame 
(and associated 'inertial time'). 

Even in the absence of a central gravitational 
field an oscillatory gravitational field will 'disen- 
gage' a spacecraft's constant speed inertial frame, 
and its geodesic trajectory will acquire an addi- 
tional longitudinal (speed) fluctuation — relative 
to a some systemic inertial reference frame. 

7.2 Agenda 

The aim of Section is to show that the (short 
term 33 ) through-time Pioneer anomaly [a p (£)] can 
alternatively be seen as resulting from a Fourier- 
like summation and superposition of 'sinusoidal' 
waves 34 , i.e. gravitational field undulations (of am- 
plitude Aa) upon the gravitational field. The re- 
mainder of Section seeks to show, by way of total 
energy concerns, that the average Pioneer acceler- 
ation (op = fflp(t)) (over a long time) is described 
by: 

a 2 P = £(Aa) 2 or a P = ^S(Aa) 2 

whereas the through-time the instantaneous Pio- 
neer acceleration [<ip(t)] may be expressed by way 
of: 

[a p (i)] 2 = £ [2 Aa 2 cos 2 (cut + (/>)] or more formally: 

n 

«pW = {£[ 2 Aof cos 2 (wit + ^)]}s 

i=l 

where <f>i is the initial (or epoch) phase of a partic- 
ular moon's wave, with orbital angular velocity Ui 
and wave amplitude Aa^ also moon specific. 

The relationship between the energy of a spher- 
ical acceleration wave upon the gravitational field 
affecting a 'point' mass through time, and the os- 
cillatory response of the moving mass to the wave 
shall be our prime concern. Recall that (systemi- 
cally) A "a* and Alt are 180° or tt c out of phase. 

32 The term "systemic" is preferred to "systematic". 
33 Especially 1 day, 5 day, or 10 day averages. 
34 Recall that 'sinusoidal' means 'harmonic' i.e. some func- 
tion built up from sinusoidal and cosinusoidal waveforms. 
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7.3 Rayleigh's energy theorem or 
Parseval's theorem 

Rayleigh's energy theorem, sometimes also known 
as Parseval's theorem, is an energy-conservation 
theorem. It says that a signal 35 contains the same 
amount of energy regardless of whether that en- 
ergy is computed in the space/time domain or in 
the Fourier/frequency domain. For a continuous 
function it is: 

/ 2 (t)dt=/ lg(f)l 2 df (1) 

-oo J —oo 

Thus, the sum (or integral) of the square of a func- 
tion is equal to the sum (or integral) of the square 
of its Fourier transform 36 . Bracewell [201 p. 120] 
notes that: "... the theorem is true if [only] one of 
the integrals exists." 

Rayleigh's theorem, by way of "the power theo- 
rem" , is also applicable to the rate of energy trans- 
fer, or power, of a 'signal' [201 p. 121]. A power 
based application follows in Section ITU The rela- 
tionship between frequency line width and Fourier 
transform amplitude in the energy Rayleigh iden- 
tity, gives way to a clear distinction between fre- 
quency and energy in the power form of the identity 
that follows. 

By way of the Rayleigh identity, the signal co- 
exists in two domains: a time domain, and a fre- 
quency domain. In the application to follow, the 
signal shall also have two physical 'faces': a field 
face, and a moving (light or low mass) body face. 

7.4 Gravitational acceleration waves 
and spacecraft kinetic energy 
(single cycle, singular wave case) 

We make use of the Fourier wave based Rayleigh 
'power' theorem identity (equation ^ of Sec- 
tion EU, but apply its frequency and particle as- 
pect, to the new circumstance of spacecraft (S/C) 
geodesic motion relative to a 'systemic' (i.e. whole 
solar system) inertial frame centered at the solar 
system barycenter. The S/C's geodesic motion 'un- 
dergoes' a pure sinusoidal (or cosinusoidal) change 

35 The terminology describing Fourier transforms is biased 
towards electrical engineering interests. 

36 Based upon: http://research.opt.indiana.edu/Library/Fou 
rierBook/ch03.html, other internet sources, and Ref. |20l 
pp. 119-122]. 



in (gravitational) acceleration 37 and this coexists 
with both: a sinusoidal change in speed (Au) rela- 
tive to an equilibrium mean value and, as we shall 
soon see, a change (per wave cycle) in barycentric 
reference frame speed (Sv). 

Since Einstein's General Relativity did not suc- 
ceed in making acceleration relative, this 'real' sys- 
temic acceleration is logically capable of represen- 
tation in a (solar-)systemic inertial reference frame. 

For a single cycle of a singular (or isolated) phys- 
ical wave upon the gravitational field of: finite ex- 
tent, fixed frequency (F), and period (At), it fol- 
lows from Rayleigh's theorem that: 

/ A Vwdt= r^pdf 

Jo Jo 

if f(t) is made an even (cosine) function. Note that 
F = (At) -1 . Regarding the right hand side of the 
equation, the single cycle's energy and power are 
easily expressed in MKS units 38 . Thus, |g(/)| 2 is 
an energy with units of [Joules per cycle] and F has 
units of [cycles per second]. Noting [<?(/) | 2 is fixed 
(per full cycle), and letting /(t) = Aacos(wt), we 
have: 

/ [Aa 2 cos 2 (wt)]dt = |g(F)| 2 / df 
Jo Jo 

Bracewell j20| reassures us that Fourier transforms 
of physically real, yet finite sinusoidal waves, do 
mathematically exist. With the tiny physical un- 
dulations of the gravitational field (Aa) physically 
affecting a moving mass, we tentatively let the 
(Fourier transform) amplitude-squared term equal 
a specific kinetic energy so that: Ig^)! 2 = ^(5vf) 2 
(a physical specific kinetic energy) . Now, since the 
integral of a squared sine or cosine function over one 
period, is half the amplitude squared multiplied by 
the period; upon integration and at a given fre- 
quency, we obtain: 

l -t\a 2 L\t = l -(5v F ) 2 F 

For convenience, the F subscript on 5v is dropped, 
and we let F — > /, with / now signifying a given 

37 This acceleration is so tiny that it would be difficult to 
measure, even with a very sensitive accelerometer. 

38 Mass, kilogram, second units. The frequency considered 
is for a single cycle, but the use of MKS units requires a 
frequency with units of [cycles per second] . 
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fixed wave frequency. Noting that ^(Sv) 2 is written 
as ^Sv 2 , the identity becomes: 

l -^a 2 M= l -8v 2 f (2) 

7.5 Discussing a single (or unit) cy- 
cle of a singular (isolated) wave 

Recalling the dimensions of the identity are 
[L 2 / T 3 ] which indicates a rate of specific energy 
transfer. The term \5v 2 is thus indicative of a spe- 
cific kinetic energy transfer associated with the si- 
nusoidal oscillation for a single cycle. 

This specific, or mass independent, energy trans- 
fer rate is also being expressed as an unorthodox 
integral of squared wave acceleration (i.e. squared 
field undulation) over time affecting a body. Note 
that, at the moment, we are not directly dealing 
with an equality involving wave amplitudes of two 
physical quantities (Aa and Aw), but an equal- 
ity of specific energy transfer rates at, and of, a 
body. Note that from the spacecraft's local refer- 
ence frame it would be the barycenter that oscil- 
lated relative to the spacecraft. Additionally, note 
that special relativistic effects are negligible in the 
solar system, while general relativistic effects are 
already incorporated in the comprehensive broader 
analysis that establishes the Pioneer anomaly. 

The two domains involved for this application 
of Rayleigh's (Power) theorem apply to: firstly (in 
the time domain), wave energy transfer rate (for 
one cycle) affecting a light body in the field; and 
secondly (in the frequency domain), [to] the rate 
of oscillatory kinetic energy expressed (in a cycle) 
by the body. The former is an effective (specific) 
energy transfer rate, and the latter a (specific) en- 
ergy 'expression' rate. Both are 'power quantities' 
determined relative to a standard systemic inertial 
frame (centered) at the solar system barycenter. 

Just what is being physically represented, i.e. a 
field energy transfer rate being equal to a kinetic 
energy transfer rate is yet to be fully clarified. 

7.6 Investigating potential energy 

With the constant phase of the undulations upon 
the gravitational field propagating at the speed 
of light, the distance covered by the accelera- 
tion wave's (constant) phase in a single period of 



Jupiter's moon Io (for example) with an orbital pe- 
riod of ~ 1.77 days, is about 300 AU. With the 
Pioneer spacecraft travelling at about 2.5 AU per 
year, spacecraft are effectively stationary as far as 
the wave's velocity is concerned. Thus, the varia- 
tion in gravitational field strength may be idealized 
to simply involve the temporal cyclical variation in 
gravitational field strength, i.e. we assume changes 
in S/C distance from the barycenter are negligible. 
We further idealize the situation by assuming con- 
stant spacecraft mass. A varying (specific) poten- 
tial energy is then simply determined by way of an 
oscillating gravitational field strength at the S/C 
or light body. This is a non-central force (or non- 
Schwarzschild geometry) based potential energy. 

Unlike a central field where (specific) potential 
energy is location dependent and proportional to 
gravitational field strength, the specific energy of 
an acceleration wave is dominantly period depen- 
dent (i.e. time dependent) and proportional to the 
square of its amplitude. Thus, for the wave's energy 
over a single oscillation, from equation [21 it follows 
that: 

— — A 2 Ay- 2 

& wave — 2 

This may be thought of as the specific undulation 
energy per cycle of the ('surface' undulatory) grav- 
itational field influencing a moving body. 

At an essentially fixed barycenter distance, in ad- 
dition to the location based (primarily) static (or 
non-oscillatory) potential energy of a body, we now 
also recognize a time based oscillatory specific po- 
tential energy associated with a wave field undula- 
tion — upon the Sun dominated gravitational field 
of the solar system. 

7.7 Relating wave potential energy 
to kinetic energy (single wave, 
singular cycle case) 

We now conceptually examine how this additional 
energy carried by acceleration waves upon the grav- 
itational field alters the expression for kinetic en- 
ergy of a body in motion. We seek to confirm the 
specific kinetic energy associated with a single os- 
cillation is \5v 2 . 

Recalling Aa = loAv we observe how the rel- 
ative oscillatory motion (Ai>, which is also some 
function of 8v, see Section l7.8|) . and potential en- 
ergy fluctuations (some function of Aa) are 180° or 
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7r radians out of phase. Just how the expression for 
specific K.E. per cycle ^(Sv) 2 , which is a response 
to specific P.E. fluctuations upon the field, is to be 
understood is now addressed. 

Total energy considerations are different with 
type-2 undulations upon the gravitational field. 
When a sinusoidal acceleration wave causes a mov- 
ing mass to display a sinusoidal (output) motion, 
it is necessary that the inertial mass 'carrying ca- 
pacity' associated with the wave's specific potential 
energy, is just sufficient to 'excite' the moving ob- 
ject, of given mass, into its sympathetic oscillatory 
motion response. Thus for the spacecraft: 

(P.E. wavc ) s/c = im s/c Aa 2 At 2 

If this is the situation, the two (inertial) masses 
are seen as common and equal for the wave's effect 
upon the moving mass. Hence, we may discuss ei- 
ther: specific energy equivalence (as in equation [5J 
or simply energy equivalence; with the two modes 
of equality not being different when resonance is 
occurring below the cut-off mass (m c ) — of a par- 
ticular wave at a particular location. Note that be- 
low the cut-off mass (m c ), the effect is independent 
of mass. The cut-off mass itself displays a different 
type of mass dependence, so that this m c is a func- 
tion of the wave's origin and evolution in space 39 . 

Relative to the inertial barycentric reference 
frame, 'light' bodies (e.g. spacecraft) are in both: 
translational motion and a cyclic (line-of-sight) mo- 
tion. Since the local frames are not mechanically- 
forced, unlike the simple harmonic motion of a 
spring, there is a specific force (i.e. an acceler- 
ation) restoring back to, and driving away from, 
the equilibrium position. A body's kinetic en- 
ergy (K.E.) 'reservoir' must be drawn upon to ap- 
pease the field's oscillatory demands. It is hypoth- 
esized that: the K.E. given to the Pioneer space- 
craft, post-planetary encounter, is being partially 
directed by P.E. oscillations (each and every cycle) 
into a 'non-productive' oscillatory response motion 
around an equilibrium value 40 . In terms of transla- 

39 The model employs E wave = im c Aa 2 At 2 at an ini- 
tial reference radius, to represent the wave's total energy. 
Ewave equals the 'gravito-quantum' excess energy (E eX cess) 
'released' to the gravitational field for (some) third bodies in 
systemic three-body motion — under the influence of weak 
gravitational fields. 

40 The violation of conservation of energy, by way of an 



tional kinetic energy, this oscillatory expression of 
K.E. is lost at the rate of ^m s / c (Sv) 2 per cycle. 

7.8 Relating motion shortfall to the 
amplitude of sinusoidal speed 

Returning to the equation (J2J equality at the end of 
Section ITU and substituting now for: Aa = ojAv 
(from Section [6.11 . At = 1// and / = oj/2tt we 
obtain: 

I( wA «) 2 A = W(f ) 

Z Ul Z ZTT 

which simplifies to give Av = Sv/ 2tt or: 

Sv = 2t:Av (3) 

Thus, the total loss of (solar-) 'systemic' steady 41 
speed in a single cycle (Sv), relative to predicted 
steady speed, equals 2tt times the sinusoidal ampli- 
tude of speed variation (Av). Note that this sinu- 
soidal amplitude in speed is about a non-constant 
mean value because kinetic energy redistribution 
acts to slow the spacecraft's steady-translational 
speed; so that at the beginning and end of the cycle 
the two speeds, relative to the (essentially fixed) 
barycenter, are different (w nna i — initial = — Sv 
where Sv > 0). Finally, note that \Slf\ — Sv, and 
that we treat Slf predominantly as simply a scalar 
loss of speed. 

7.9 Steady and unsteady energy of 
motion (singular wave case) 

The total kinetic energy, possessed by the Pioneer 
spacecraft, can be seen as equal to a steady kinetic 
energy component plus an unsteady kinetic energy 
component. Note that for multiple waves from mul- 
tiple moons, of distinct period and amplitude, there 
will be a number of coexisting and superpositioned 
unsteady component terms. 

Importantly, one needs to appreciate that in the 
absence of any sinusoidal speed variations all the 

asymmetrical quantum indeterminacy collectively shared by 
each and every atom in a moon, expressed in the guise of 
a field undulation, has eventuated in this further (or sec- 
ondary) new phenomenon. 

41 'Steady' here means non-oscillatory as compared to os- 
cillatory (i.e. unsteady). Naturally, the steady motion may 
be: linear (radial), rotational (i.e. circular or elliptical), hy- 
perbolic or parabolic. 
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kinetic energy of a moving body is steady (i.e. non- 
oscillatory) . This being the case with general rel- 
ativistic orbital mechanics, i.e. how the solar sys- 
tem's gravitation is currently modelled. For the 
Pioneer anomaly (in the far outer solar system) 
small non-conservative effects, such as radiation 
pressure and antenna reaction forces, may be safely 
neglected. 

The motion of both Pioneer spacecraft are pre- 
dominately steady, and idealized to be radially di- 
rected 42 . A sinusoid in S/C speed (of amplitude 
Av), driven by a single sinusoid in gravitational 
field strength (of amplitude Aa), will redistribute 
some of the total (kinetic) energy of a body from 
fully steady motion into an additional unsteady (or 
spectral) component of specific K.E (\5v 2 per cy- 
cle). This longitudinal oscillation around an 'av- 
erage' velocity requires kinetic energy, and this 
requirement results in a shortfall in translational 
speed of Sv per oscillation — relative to the pre- 
dicted condition of fully 'steady' motion. This addi- 
tional component of K.E. (longitudinal oscillatory 
for the Pioneer spacecraft) is 'expressed' each and 
every cycle. Over a long time these Sv increments 
will sum to significant levels. Over a single cycle 
of a singular (effectively monochromatic) wave we 
have: 

(AK.E.) stca dy = (K.E.) uns toady = ^m s/c (Sv) 2 

The P.E. oscillations act to continuously 'erode' 
the spacecraft's steady kinetic energy, thus caus- 
ing a slowing as compared to predicted motion. 
This slowing effectively mimics a non-conservative 
force-like drag effect upon (the hitherto assumed) 
essentially frictionless motion. Think of it as an 
acceleration- wave based drag if you like, that 'chips' 
away a fixed amount of steady K.E. each and every 
cycle, to 'partner' the P.E. oscillation. Hence, total 
kinetic energy, i.e. steady plus unsteady K.E., is 
also chipped away. Remember this effect only ap- 
plies to 'light' bodies below the distance dependent 
cut-off mass (m c ) of each specific wave. 

7.10 Interim Remarks 

Consequently, the identity between Aa and Sv 
(equation 0) derived from Parseval's theorem 

42 Their outer solar system hyperbolic motion is idealized 
to have an infinite radius of curvature. 



(equation appears to have a justifiable physi- 
cal application, i.e. conceivably explaining the Pi- 
oneer anomaly. New physics, beyond the notion of 
the problematic extra force hypothesis 43 , is arising 
out of a contemporary interpretation of established 
mathematics and physics. It is the 'mechanism' 
behind the establishment, and the associated quan- 
tification of the wave energies and m c magnitudes, 
that remains in need of rigorous scientifically expli- 
cation. 

In this write up the waves' existence is simply 
being hypothesized and then scrutinized. Partial 
aspects of a more detailed theoretical model have 
been included only to facilitate this process. 

7.11 Multiple acceleration waves 
upon the gravitational field 

We are now in a position to examine the effect of 
the superposition of acceleration waves, emanating 
from a number of moons, upon a spacecraft. Let us 
continue to idealize the situation by letting the S/C 
be in pure radial motion at a position far outside 
the solar system, and by letting the total observa- 
tion time (i.e. the total data interval) be very long. 

Note that the physical mechanism hypothesized, 
over say 100 years, has effectively fixed values of: 
Aa, Aw, Sv and At = 1// for any given moon. 

From equation (i.e. Sv — 2ttAv) and u> — 2irf 
it follows that for a singular wave: 

(Av)u = (Sv)f 

where (Sv)f, or if you prefer (Sv • /), is the rate of 
speed shortfall (in MKS units) for a singular wave, 
over a single wavelength. Now since Aa = luAv. 

£(Aa) 2 = S(An-w) 2 
= Z(Svf) 2 
= ^(Sv/At) 2 

Since we are dealing with energy transfer, we must 
continue to work with squared quantities, in this 

43 An additional force, i.e. additional spacetime curvature, 
cannot explain the Pioneer anomaly's apparent violation of 
the principle of equivalence, since the orbits of planets over 
a considerable period of time would indicate the presence of 
such a 'generally applied' force (i.e. the force influences all 
masses). A similar situation applies to the orbits of long- 
period (heavy-body) comets 1211 . Additionally, large outer 
solar system asteroids (of diameter >~ 3 km) will not ex- 
hibit the Pioneer anomaly 1221 . See Section l8.2l 
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case Aa. Thus, ^S(Aa) 2 or [S(Aa) 2 ]3 is used to 
represent the overall sum of all acceleration ampli- 
tudes. This sum equals long-term average (addi- 
tional) acceleration opor a p (t). Thus, in the ideal- 
ized circumstances discussed, it is feasible to write: 

a P = ^(t) = \/£(Aa) 2 = [£(Aa) 2 ] i 

Now since Doppler tracking measurements, and 
a raft of supporting science, give a value of anoma- 
lous speed shortfall (AVp), that over a medium- 
long period of time (Tj,) indicates the action of 
an apparently constant anomalous acceleration, we 
may establish that: 

lim [S(Aa) 2 ]i = lim [Z(5v/At) 2 ]i 

t^largc t— >largc 

AVp 

= dp 

= constant 

Over extended periods of time (i.e. 11.5 years) we 
have (on average) a 'headline' constant result so 
that: ap — £(Aa) 2 = constant. Whereas in the 
short term the a p (t) measured has a stochastic-like 
variation 44 about the long term fixed mean value 
(dp). Thus, upon (hypothetical^) removing all sig- 
nal noise, it is proposed that: 

[a. p (t)] 2 = £[2Aa 2 cos 2 (wt + <j>)\ 

best matches the observational data to the hypoth- 
esized physical situation 45 . Note that the posi- 
tional variation in the additional (Pioneer) accel- 
eration [a p (x,t)] has not been included, as yet, in 
this idealized model; where x is distance from the 
barycenter (i.e. from the systemic inertial frame's 
'origin'). 

44 Actually, the real lunar-based variation is determinis- 
tic if the wave amplitudes and cut-off masses are well de- 
scribed, and the planetary positions are accurately known. 
Even though the acceleration waves are effectively fixed in 
amplitude and frequency (over long-terms), the motion of 
their host planets and the motion of a spacecraft (or 'light' 
body), will alter the effect at the spacecraft through time. 
Thus, the direction-based variations in amplitude produce 
a deterministic time series that coexists with the effects of 
stochastic Doppler signal noise. 

45 The factor of 2 is necessary, (tentatively) in order to 
make the through-time changes of the sum of the solitary 
waves' individual contributions of 5vi(t) — equal to [a p (t)]. 



For a radially (outward) directed spacecraft, at a 
given position and a given time, we may schemati- 
cally express the total (inward) acceleration acting 
on a 'light' body as: 

<hotai(%, t) = ogr(i, t) + Op(x, t) - a rad (x, t) + ... 

where a rac i is (specific) outward time dependent so- 
lar radiation force. Note that both the: total in- 
stantaneous acceleration and this additional (or Pi- 
oneer anomalous) acceleration are now functions of 
both location and time. Additionally, in a three 
(or more) body system ogr, or standard gravita- 
tional thcorization, is also time dependent. The 
non-determinism of these circumstances produces 
a need for ephemerides. 

From the above it is also clear that the addi- 
tional (or Pioneer anomalous) acceleration (a p ) is 
independent of overall spacecraft speed. 

7.12 Discussion, further remarks, 
and relaxing the idealization 

This section has sought to illustrate that the Pio- 
neer anomalous acceleration may (alternatively) be 
hypothesized as resulting from the effect of a su- 
perposition of 'sinusoidal' wave undulations upon 
a pre-existing gravitational field. Observations and 
the author's wider model imply the undulations are 
lunar based, with only some (5 out of the 7 major) 
moons 46 contributing for kinematical and geomet- 
ric reasons. Geometry is involved by way of both 
spacetime curvature and quantum mechanical geo- 
metric phase. 

The long term mean amplitude of all the waves 
together is fixed (or constant) but over shorter time 
periods there is necessarily variation around the 
long-term mean acceleration. The statistical vari- 
ance of observations of a p (x,t) has been greatly 
moderated (i.e reduced) by the orbital resonances 
of the Galilean moons of Jupiter. 

When spacecraft are in the 'mid' (10 AU) to 
outer solar system the spherical wavefronts are not 
approximately orthogonal to the direction of S/C 
propagation, and directional cosines play a signifi- 
cant role 47 . The anomaly will vary with the posi- 

46 For a variety of reasons the other prograde spin-orbit 
coupled moons of the solar system may be neglected. The 
main reason is that: Aa is proportional to total lunar mass. 

47 Within the orbit of Saturn things are different again, 
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tions of Jupiter and Saturn relative to the space- 
craft. The effect goes to a slight asymptotic max- 
imum as a spacecraft leaves the solar system and 
its heliocentric radius gets very large 48 (so that di- 
rection cosines between the direction of the moving 
wavefronts and a body's trajectory approach one, 
i.e. unity). Thus, where x is distance from the 
barycenter: 



ap = a p (x,t) ~ constant if a; — » large 

Large x is > 20 AU. Surprisingly, due to direc- 
tion cosines and the possibility of obtuse angles 
inwards of Saturn, we have in general (especially 
for x < 10 AU): 

ap =/= constant for all x 

which is supported by early Pioneer 11 observa- 
tions Fig.7, p. 18]. 

A second look at the anomaly in the Galileo 
spacecraft's cross solar system journey out to 
Jupiter (Dec 1992- July 1995) might clarify this is- 
sue. Note that a trajectory or path-based hypothe- 
sis, as compared to a Sun-based effect, involves dis- 
tinctly different angles, and hence different anoma- 
lous acceleration magnitudes, at these low AU val- 
ues — of 1 to 5 astronomical units. 

The author's fuller model indicates that at all 
times the spacecraft's mass is many orders of mag- 
nitude below the 'cut-off' mass at which point a 
body no longer resonates with at least one of the 
undulations upon the gravitational field 49 . This is 
touched upon in Section HT2l 



with acceleration waves (of significance) being able to re- 
tard a body's motion at obtuse angles to each other, i.e. in 
opposing directions. (Recall Section 16.51 ) 

48 The author's fuller model yields an asymptotic value of: 
a P = ^/[E(Aa) 2 ] £3 9.16 X 10" 8 cm/s 2 . Once direction 
cosines (both in the plane of the ecliptic and inclined to 
it) are incorporated, the value of a p for Pioneer 10 (1987 to 
1998) reduces to ~ 8.8 X 10 — 8 cm/s 2 . This is embarrassingly, 
yet significantly, close to the stated result of (8.74 ± 1.33) X 
10 -8 cm/s 2 given by Anderson et. al. 

49 The existence and magnitude of the cut-off masses, as a 
function of distance, indirectly has its origin or source in the 
mass dimension associated with Planck's constant, and ex- 
pressions for maximum indeterminacy in quantum mechani- 
cal situations. (Without a fuller model of acceleration wave 
'generation', this qualifying footnote is insufficiently concep- 
tually supported and therefore it may be ambiguous.) 



8 Further applications of the 
hypothesis 

The hypothesized existence of: lunar generated, 
constant amplitude, first order fluctuations upon 
the gravitational field also potentially affects, in a 
beneficial way, other anomalous circumstances in 
our solar system. 

8.1 The Earth fly by anomaly 

Assuming the reality of this anomaly, the preceding 
hypothesis has a promising, easy and direct quali- 
tative application to the anomalous increase in the 
velocity of spacecraft associated with some, but not 
all, Earth fiybys [23! El 

The kinetic energy of the inbound spacecraft will 
have tiny additional (longitudinal) oscillatory com- 
ponents of kinetic energy that are overlooked 50 . At 
Earth gravity assist encounter, in the planet's ref- 
erence frame, total S/C energy remains constant 
(i.e. is conserved) P5, p. 449]. Then, post encounter 
either: the change in trajectory orientation rela- 
tive to Jupiter and Saturn reduces the total of the 
oscillatory components, or alternatively (and less 
likely), any longitudinal oscillatory component is 
yet to fully reestablish itself. Subsequently, an ex- 
cess (steady) kinetic energy over predictions will be 
evident on occasions, and the effect will vary de- 
pending upon the pre- and post-encounter trajec- 
tories relative to the positions of (the moons hosted 
by) Jupiter and Saturn. 

It is very difficult for any model to explain an 
anomalous increase in kinetic energy Other than 
to cite observational errors, the hypothesis pre- 
sented herein is conceivably the only reasonable 
non-systematic physical explanation possible. 

8.2 Early solar system history 

There are at least two problems concerning the 
outer solar system's very early history 51 for which 
a new slowing or 'braking' mechanism, in regions of 
very weak gravitational fields 52 , may prove benefi- 

50 Numerous pre-encounter trajectory manoeuvres may 
possibly interrupt and/or overwhelm this oscillatory com- 
ponent of kinetic energy. 

51 After lunar formation and spin-orbit coupling had been 
attained, around their respective 'host' planets. 

52 Note that this excludes planetary ring systems. 
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cial. Like the (predominantly) radial Pioneer case, 
orbital motions (in the outer solar system) are sub- 
ject to a predominantly radially directed oscilla- 
tion, that in this case is now orthogonal to the tra- 
jectory of their motion. As for the Pioneer anomaly, 
this represents an unsteady expression of kinetic 
energy. The inward acceleration of the Pioneer 'ra- 
dial' motion, for an orbiting body becomes a decel- 
eration of magnitude 8.74 x 10~ 8 cm/s 2 effect. Over 
a quarter of a million years this would produce an 
approximately 7 km/s slowdown of orbiting 'light' 
bodies 53 , while leaving 'heavy' bodies (and electro- 
magnetic radiation 54 ) totally unaffected. Note that 
Neptune's average orbital speed is 5.43 km/s. 

For 'big rocks' (below the cut-off mass transi- 
tion), this new physical 'method' of braking light 
bodies and hence altering their motion, is much 
'quicker' than any conceivable braking by existing 
conventional drag effects. 

The migrating planets hypothesis that is invoked 
to explain the too rapid formation of the ice giants 
(Uranus and Neptune), in their present locations, 
may also possibly be overcome. This hypothesis is 
not without its concerns 55 and criticisms [2U |27|. 
If all small primordial bodies, in very weak grav- 
itational fields, underwent 'wave-braking' by the 
aforementioned mechanism then they would have 
spiralled into outer planets or the inner solar sys- 
tem. Compared to other 'braking' mechanisms this 
happens very quickly and very early in the solar 
system's history. Naturally, this mechanism con- 
tinues to exert its influence today, with collisions 
of bodies generating smaller bodies 'continuously'. 

There is also an apparent total lack of 'small' 
comets observed in the solar system less than about 
lkm in diameter [281 Fig.2], or ~ 5.2 x 10 n kg (as- 
suming a comet density of 1 x 10 3 kg/m 3 ). This 
implies that primordial small comets were not flung 
out, from the inner solar system, to the Oort cloud 
long-ago for some reason. The author's fuller model 

53 For bodies inward of Uranus the deceleration will be a 
bit less, or slower, than these values. 

54 The speed of massless electromagnetic radiation pho- 
tons also remains (necessarily) unaffected. The dual-aspect 
acceleration-waves on the gravitational field, involving both 
an oscillatory acceleration magnitude component and a 
scalar inertial mass component, are restricted to only in- 
fluencing (light) bodies of non-zero mass. 

55 For example, distant Neptune's orbit, although rela- 
tively easily disturbed, is 2nd only to Venus in its circularity 
(i.e. low eccentricity). 



produces different cut-off masses (m c ) for each of 
the five moons previously mentioned — that are 
considered to dominate the Pioneer anomaly. At 
2.5 AU from the lunar sources, these m c values 
range between 3.4 x 10 12 kg (Titan) to 4.5 x 10 13 kg 
(Callisto). Whereas at 20 AU from the lunar 
sources, these values go from 6.6 x 10 9 kg (Titan) to 
8.9 x 10 10 kg (Callisto). Observational evidence 
Figs.l and 2] implies a depletion in comets, as com- 
pared to expectations, beginning at about 5 km di- 
ameter, or a mass of ~ 6.5 x 10 13 kg. Bearing in 
mind that primordial comets (orbiting in the plane 
of the ecliptic), exist over a range of heliocentric 
distances with somewhat random orbital motions; 
the cut-off masses derived appear to qualitatively 
'match' this situation in a very promising way 56 . 

In a similar vein, Kuzmitcheva and Ivanov |28j 
conclude: "The lack of small craters on [the as- 
teroid] Eros is the first observational evidence of 
a possible paucity of smaller bodies in the Main 
Belt ." Additionally, there apparently exists an 
'abrupt' edge to the present (and presumably the 
long-past) solar system by way of an abrupt outer 
edge to the Kuiper belt. This was the exact oppo- 
site of what planetary scientists had envisaged 30 . 

Finally, the evidence for this ongoing braking ef- 
fect on 'light' bodies (of non-zero mass) appears to 
be supported by the contemporary representation 
of a kink in the linear relationship between the cu- 
mulative number of Earth crossing asteroids vs. the 
diameter of these bodies (see the log-log diagram: 
Fig. 3 of Ref. 28 ) . There is one straight line for 
big primordial bodies, and the other for 'fall-out' 
from collisions of these bigger bodies. The kink 
joining the two possibly illustrates both: the size 
level at which bigger primordial bodies start to be- 
come sparse, and (at the other end of the kink) the 
upper size limit of impermanent collision debris. 

9 Concluding discussion 

Observational evidence of the 'through time' be- 
haviour of the Pioneer anomaly a p {t) is open to a 
new interpretation by way of hypothesizing a num- 
ber of tiny (sinusoidal) undulations (Aa) on the 

56 The author's model has not been based upon this evi- 
dence in any way. The observations were found to support 
the model after it was produced to possibly account for the 
Pioneer anomaly's observational evidence. 
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gravitational field; each with its own fixed ampli- 
tude and frequency. Ignoring direction cosines: 



Op(t) = 



^][2Aa 2 cos 2 (^< + &)] 



and, on average over a long period of time, e.g. 
30 years (roughly the orbital period of Saturn): 



ap = a p (t) = \/E(Aa) 2 

The acceleration-wave hypothesis is supported, in 
this write up, by partial aspects of a new model that 
seeks to explain the observational evidence imply- 
ing at least three anomalous behaviours in our solar 
system. The hypothesis, at this stage, is not easy to 
accept, because the new physical model discussed 
has not been fully established. This needs to be 
presented to support and compliment the predomi- 
nantly general presentation of a prospective model 
given in this article. The scope of the hypothesis is 
unlike other 'singular solutions' offered to explain 
a 'constant' Pioneer anomaly. Importantly, in this 
new approach the apparent violation of the Equiv- 
alence Principle is accepted and incorporated at a 
fundamental level in the model. 

The Pioneer (acceleration) anomaly essentially 
becomes an oversight in predicted spacecraft mo- 
tion, by way of omitting the effect of certain 
wave-like undulations in the gravitational field's 
strength. A shortfall in motion results which is 
effectively a new non-conservative (specific force) 
'deceleration-wave- drag' effect. The most impres- 
sive evidence for this effect, that is restricted to 
bodies whose mass is non-zero and below a certain 
cut-off mass zone 57 , is the Pioneer anomaly. 

This article has sought to show how an ongoing 
braking effect, applicable only to light bodies has 
sculptured the contents of our solar system from its 
very earliest stages. Somewhat ironically, it may be 
conjectured that: were it not for the acceleration- 
wave-braking mechanism proposed herein, the Pi- 
oneer (and other) spacecraft may not even have 
passed safely through the asteroid belt on their 
journeys to Jupiter and beyond, many years ago, 
because it would remain strewn with collision frag- 
ments of all sizes. 
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